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Abstract 

We report on a macroscopic version of the single-electron transistor (SET), which we call the 
soliton tunneling transistor (STT). The STT consists of a gate capacitor coupled to a NbSea crystal 
with a charge density wave (CDW). We find that the current-voltage characteristic of an STT is 
periodically modulated by the gate voltage, as in the SET, except that the measured periodicity 
corresponds to a macroscopic displacement charge. These results appear to be consistent with time- 
correlated quantum nucleation of solitons and antisolitons [see Phys. Rev. Lett. 84, 1555 (2000)]. 
We discuss how the microscopic degrees of freedom within the condensate might enable quantum 
behavior at high temperatures, and report on preliminary modeling studies using a coupled-phase 
Hamiltonian to interpret our results. 
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I. BACKGROUND AND MOTIVATION 



A wide class of nonperturbative phenomena in field theory can be understood in terms of 
quantum tunneling. A well-known example is the quantum decay of the false vacuum,§ which 
has been of broad scientific interest in cosmologyi'i and other fieldsi for over two decades. 
In three dimensions, the boundary between the bubble of true vacuum and the surrounding 
false vacuum is a type of topological defect known as a domain wall. A variety of topological 
defects in condensed matter systems have been proposed to nucleate via quantum or thermal 
fluctuations. These include vortex-antivortex pairs and vortex rings in superconductorsj'i 
superfiuids,@'i and Bose-Einstein condensates;! dislocation pairs in Wigner crystal^Sl and 
vortex lattices;tiJ phase-slip vortex rings in charge density waves (CDWs);00 charge (or flux) 
soliton-antisoliton pairs in density wavesEj'Hj (or Josephson junctionslii); and soliton- 
like domain walls surrounding cigar-shaped bubbles of true vacuum in three-dimensional 
CDWs.lli The quantum decay of the metastable false vacuum of a scalar field 0, accompanied 
by the creation of solitons and antisolitons in (1+1) dimensionJUi (or soliton domain 
walls in (3+1) dimension^), has been studied extensively in quantum field theory. 

A CDW is a condensateil in which the electronic charge density in a quasi- 1-D metal is 
modulated, p{x, t) = po{x, t) + pi cos[2A;i?x — 0(x, if:)]. Here po{x, t) contains the background 
charge of the condensed electrons, and an excess or deficiency of charge proportional to 
±90/9x. Weakly pinned density waves (DWs) and weakly-coupled Josephson junctions 
(JJs) can be approximately described by a phase, 0(x, t), in a sine-Gordon (s-G) potential. 
They are dual in that the roles of charge and flux are interchanged, as well as those of current 
and voltage. The current in a density wave is / = {Qo/2tt) d(f)/dt, where Qq ~ 2Nche and 
Nch is the number of parallel chains, whereas the voltage across a JJ is = ($o/2vr) dcp/dt, 
where $o = h/2e. Charge (flux) solitons in a density wave (J J) carry a charge (flux) 
of ±Qo{±^q) (see Table I). The width of a Josephson vortex is roughly the Josephson 
penetration length, Aj oc J~^^'^ whereas, in a DW, the soliton width is Aq = cq/cuq, where cq 
is the phason velocity and ujq is the pinning frequency. Thus, Aq will increase with decreasing 
impurity concentration (as Uq decreases), and may approach the distance between contacts in 
extremely pure samples. This is equivalent to approaching the short junction limit, L < Aj, 
in a JJ, where the V-I curves become significantly less rounded. 
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due to the electrostatic energy of the charged sohtons and antisohtons, for the creation 
of sohton-antisohton {S — S') pairs in density waves. A more recent paperll proposes an 
analogy to time-correlated single electron tunneling to explain the observed lack of DW 
polarization below threshold, and to interpret key features of DW dynamics, such as coher- 
ent oscillations and narrow-band noise, above threshold. Observations of Aharonov-Bohm 
(A-B) oscillationsii in the magneto-conductance of CDWs in NbSes crystals with colum- 
nar defects provide compelling evidence for quantum transport. Additional evidence for 
quantum behavior is found in rf experiment sH'il'ii that show good agreement with photon- 
assisted tunneling theory, and the observed Zener-like I-V characteristics.^ This Zener-like 
behavior is expected in the quantum picture when L ^ Aq and the I-V curves are rounded. 
The linear I-V curves seen in the extremely pure NbSes samples suggest a possible Coulomb 
blockade mechanism, the dual of which is also suggested by the by the linear V-I curves 
seen in cuprate JJs in the short junction limit (L < Aj). 

Evidence against classical depinning includes bias-independent rf and microwave re- 
sponses below threshold in CDWs^'0 and Wigner crystals,!! and the observed small phase 
displacements of chargeiiS and spinil density waves below threshold. These experiments 
strongly suggest that, as the electric field is increased, S — S' ot dislocation pairs are cre- 
ated long before the classical depinning field is attained.ii Moreover, attempt^ using a 
scanning tunneling microscope (STM) to directly observe either displacement of the CDW 
below threshold or sliding above threshold have been unsuccessful. The apparent lack of 
sliding seen in STM experiments, the jerky dynamics revealed by NMR experiment s,ii and 
the mode locking observed at high drift velocities^ all suggest that the CDW spends most 
of its time in the pinned state even above threshold. 

The quantum interpretation of the threshold field, as a pair-creation threshold due to 
Coulomb blockade, is motivated by Coleman's paperS on soliton pair-creation in the massive 
Schwinger model. A pair of S and S ' domain walls with charges ± Qo produce an internal 
field of magnitude E* = Qo/eA, as shown in Fig. 1, where A is the cross-sectional area. 
When a field E is applied, the difference in electrostatic energies of a state with a pair of 
separation / and of the "vacuum" is AU = \ eAl\{E ± E*f - E'^] = Qol[^E* ± E], which is 
positive when \E\ < ^ E*. Conservation of energy thus forbids pair production for fields less 
than a quantum threshold, Et = ^E* = Qo/2eA, which appears to be about two orders of 
magnitude smaller than the classical depinning field in NbSe3.ii The observed universality 
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relation, eEt ~ eNch/A^ thus arises quite elegantly in this model.0 

A density wave between two contacts behaves as a capacitor with an enormous dielectric 
constant (Fig. 1). The initial charging energy is Q'^/2C, where Q is the displacement charge 
and C = sA/L. We define 9 = 2ttQ/Qq = 2ttE/E* = ttE/Ex and note that a displacement 
near the middle creates a non-topological kink-antikink pair with charges ±(0/27r)(5o, if 
= at the contacts. The washboard pinning and quadratic charging energies can then be 
written asS 

rL 

U[(j)]= dx {up[l- cos (j){x)] + Uc[e -(j){x)]'^} (1) 

J 



where the first and second terms represent the pinning and electrostatic charging energies, 
respectively, and where Up ^ Uc for NbSes.B If the system starts out in its ground state, 
conservation of energy will prevent tunneling when the applied field is below threshold, 
9 < -K [E < Et), as illustrated in Fig. 2. However, when 6 exceeds vr what was formerly 
the true vacuum becomes the unstable false vacuum. A bubble of true vacuum, with soliton 
domain walls at its surface, then nucleates and expands rapidly (Fig. After n solitons of 
charge Qq (and antisolitons of charge —Qo) have reached the contacts, the charging energy 
becomes 
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This series of piecewise parabolas is similar to the charging energy of a single-electron tunnel 
junction, except that Qq now represents a macroscopic charge. 

A single-electron transistor consists of a gate capacitor Cg coupled to an island 

electrode between two small capacitance tunnel junctions in series. The gate voltage modu- 
lates the I-V curves between the source and drain electrodes, with a period e in displacement 
charge, Qg = CgVg. The displacement charges Qi^2 across the two tunnel junctions are re- 
lated as Q2 = Qi + Qg + qo, where go is a phenomenological offset charge induced during 
cooling.ii The SET is related by charge-fiux duality to the dc SQUID. The critical voltage 
across an SET is a periodic function of Qg, whereas the critical current across a SQUID is 
periodically modulated (with period $0) by the fiux $. 
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The model discussed above suggests that it may be possible to demonstrate a macroscopic 
version of the SET, by attaching a gate capacitor to an island electrode near the center of 
a quasi- 1-D crystal with a density wave. The displacement charge induced by the gate 
electrode would then periodically modulate the total critical voltage between the source and 
drain electrodes. Ideally, in the absence of any shunt conductance, the periodicity of the 
gate displacement charge might be expected to be ~ Qo- However, screening by the normal, 
uncondensed electrons will tend to reduce the effectiveness of the gate which, unlike the 
source and drain contacts, cannot be driven by a current source. The displacement charges 
across the two segments of the crystal will be related as Q2 = Qi + PiQg + lo), where 
f] ~ exp{—Leff/Ls) <^ 1 and Lg ~ 2-7?/ /e^ is a screening length. The total charging energy 
of the two segments, in this idealized model, will then be 

2Ci 2C2 87r2 \ 2Ci 2C2 ] 

where Ci and C2 are the capacitances of the two segments separated by the island electrode. 
The analogy to the SET suggests that a gate voltage might modulate the I-V curves between 
source and drain contacts, with a periodicity AV^ ~ Qo/PCg- The gate capacitance Cg, and 
hence the attainable displacement charge Qg, may have been too small to observe non- 
monotonic behavior in previous experiments,^ in which a gate electrode was fabricated 
directly on the crystal to form a MOSFET-like structure. The soliton tunneling transistor 
(STT), discussed in the next section, employs a much larger, l-fiF, gate capacitor coupled 
to an NbSea crystal, and exhibits non-monotonic behavior. 



II. EXPERIMENTAL MEASUREMENTS 



Single crystals of NbSes were employed in the experiments reported here. This material 
forms two independent CDWs, at Peierls transition temperatures of 145 K and 59 K,E1 
respectively. The Peierls gap opens up over most of the Fermi surface (FS) below the lower 
transition, but leaves a small portion of the FS intact, so that a significant concentration 
(~ 6 X 10~^^ cm~^) of normal, uncondensed carriers remain down to low temperatures. The 
geometry used in our experiment is illustrated in the inset to Figure 3, where the width 
of the crystal is exaggerated for clarity. The NbSea crystal was placed onto an alumina 
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substrate with a series of evaporated, 25-/im wide gold contacts. The substrate was thermally 
anchored to a cold-finger in the vacuum shroud of an open cycle helium flow cryostat and the 
temperature was controlled using a Lake Shore temperature controller attached to a heater 
coil wrapped around the cold-flnger. A Keithley programmable dc current source injected 
the current between two contacts, which were bonded to the crystal near the ends using silver 
paint. The "source-to-drain" voltage was measured between two additional gold contacts, 
as illustrated, and a 1-fiF gate capacitor was attached to the center gold contact using silver 
paint. The spacing between contacts along the crystal was 500 fim center-to-center, and the 
gate capacitor was kept inside the cryostat. 

We found that substantially smaller gate capacitors (as well as gate capacitors with longer 
leads) were unable to induce a periodic modulation of the I-V characteristic. Moreover, 
dc I-V (rather than differential dV/dl) measurements were necessary to avoid inducing 
a displacement current through the gate capacitor. A programmable voltage source was 
coupled to the gate capacitor via a 10-kfi resistor, which limited the current flowing through 
the crystal during changes in gate voltage when the gate capacitor either partially charged 
or discharged. The cryostat was kept inside an electromagnetically shielded enclosure, and 
Vsd was measured with a nanovoltmeter. 

The measurements were primarily carried out at 35 K. Previous "fleld effect transistor" 
experimental showed the greatest modulation at around 30 K, where the threshold fleld 
is near its minimum. We also observed the largest modulation, using our geometry, at 
comparable temperatures. We attained the best temperature stability (better than ±0.01 
K) when the temperature was set to 35 K, and thus chose this temperature for most of the 
measurements reported here. 

Figure 3 shows several plots of CDW current as a function of source-to-drain voltage, 
hdw vs. Vsd, in a NbSea crystal at 35 K, for different values of gate voltage Vg. The gate 
voltage is seen to modulate the threshold voltage in the I-V curves of Fig. 3. Figure 4 
displays plots of source-to-drain voltage Vsd vs. Vg for three values of total bias current 
above threshold. The plots exhibit roughly periodic behavior, similar to that observed in 
SETs. However, in our system, the measured periodicity AVg ~ 10 V is consistent with a 
macroscopic displacement charge AQ = CgAVg ~ 6 x 10^^ e, comparable to the charge of 
the conducting electrons between the contacts. 

The behavior shown in Fig. 4 is quite extraordinary, and appears to be consistent with 
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the soliton tunneling hypothesis. The number of parallel CDW chains, Nch-, is about 10^. 
Thus, one might estimate the screening parameter as follows: (3 ~ Qo/AQ ~ 2Nche/AQ ~ 
3 X 10^^. Assuming that (3 ~ exp(— Lgjj/Ls) ~ 3 x 10^^, and taking L^ff to be about the 
distance between contacts (~ 500 /im), yields an effective screening length Lg of about 40 
/im. However, an alternative interpretation for the observed periodicity might be that all 
of the normal electrons between the contacts participate in screening out the displacement 
charge Qo of the CDW. Thus, the observation that AQ/e ~ 6 x 10^^ is roughly the number 
of conducting electrons between contacts may not be a coincidence. Further work is needed 
to better understand the effects of screening by the normal, uncondensed electrons. 

III. MODELING STUDIES 

We believe that phase-coherent, Josephson-like tunneling of many microscopic degrees 
of freedom enables quantum transport to dominate at all temperatures below the Peierls 
transition temperature, Tp. On the one hand, the observation of complete mode locking 
with an ac source in high quality crystalJ^i shows that the phase is coherent throughout 
macroscopic regions within the crystal. On the other hand, magneto-transport experiments 
on NbScs crystals with columnar defects^ demonstrate Aharonov-Bohm (AB) oscillations 
with a periodicity of h/2e, and not h/2Ne as predicted theoreticallyil, where N is the 
number of coupled chains. These apparently contradictory results can be reconciled by 
treating a density wave as a condensate containing many quantum degrees of freedom within 
a phase-coherent volume. 

This suggests writing down a Hamiltonian in terms of the phases of individual standing 
waves created by the dressed electrons, coupled to the 2k p phonons, in the CDW condensate: 

H = J2Hn + J2 Un,n'[^ " COs(0„ - (j)n')], (4) 
n n,n' 

where 0^ = |(0nt+0nj,) , Un,n' is related to the coupling between nearby chains (or transverse 
wavevectors) , and 

Hn = EfY. Kl'^ - COS(0„.)] + E'^[(l)na ' Of] + U' [I - COs(0„t " 0ni)]. (5) 
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Here 0nT(0ni) is the phase (near the mid-point between the kink and antikink) of the standing 
wave created by a delocahzed spin-up (spin-down) dressed electron in the condensate, Una — 
—id/d(j)na is the momentum operator, and Ep is the Fermi energy. The parameter [i ~ 
Mp/m is the Frohhch mass ratio, while E'^ = Ep/ Ep, E'^ = E^/Ep, and U' = U/Ep 
represent normalized pinning, charging, and coupling energies, respectively. Note that, in a 
spin density wave, the spin-up and spin-down CDWs are out-of-phase, so the last term in 
Eq. (6) would be U'[l + cos{(t)n^ — 0ni)] in an SDW. The effective charge of 2e observed in 
the AB experiments suggests that the coupling U between the two spin components is strong 
compared to the other coupling terms, i.e. U » Unn>- We thus take U' to be nonzero in 
our calculations, and consider the case where U' >> E'^ >> E'^. Alternatively, if the phases 
4>ni and (pni are taken to be identical (i. e., if [/ —>■ oo), our model can be treated as one in 
which pairs of chains are coupled. 

In the variational approach, the energy is minimized to estimate the ground-state wave- 
functions and energies for different values of 9. Because a CDW is highly dissipative, the 
system will tend to quickly relax to its lowest energy state. In the variational method, the 
energy expectation value 



E^{^\H\^) (6) 

is minimized by setting 5E — to estimate the ground state energy. 

The wavefunctions of the system xi'Pna) are approximated as superpositions of sharply 
peaked Gaussians centered at the pinning potential minima. The state of the system, in- 
corporating the two spin components, is described as a trial function with properties that 
depend on the parameters 6^ and a, 



N 



X{(l>na) ^ 6^exp[-Q;(0„^ - 27rm)^ 

m=-N 



(7) 
(8) 



The coefficients bm asymptotically approach zero as m approaches infinity. We set = 2 
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in Eq. (7) to render the problem computationally tractable. For each quantum degree of 
freedom we choose the coefficients to satisfy the following normalization condition: 

j:hl = l. (9) 

m=-2 

We have studied the problem using of a variety of parameter values in the Hamiltonian of 
Eq. (5), with qualitatively similar results. For the results illustrated here, we take /i = 354§, 
Ep/Ep = 10'^ , Ec/Ep = 10"^, and U /Ep = 5 x 10"^. Fig. 5 shows a plot of the minimized 
energy E{6). The divergence for \6\ > 47r is an artifact of our having limiting the number 
of coefficients bm to 5, i.e. taking = 2, in Eq. (7). An exact calculation would yield a 
plot in which E would be a periodic function of 6. In Fig. 5, the energy is reduced slightly 
at the crossing points, 6 = mr, as compared to piecewise parabolas, but Eq. (2) provides a 
reasonable approximation to E{6): 

E{d) ^{9- 2-Knf. (10) 

As in the SET, the voltage across each segment of the STT is related to the charging 
energy as Vi^2 = dE/dQi^2 oc dE/d6i^2- If we use the approximation given by Eq. (10), this 
yields a sawtooth function, which can be expanded as a Fourier series: 

Y °° f— ll" 

V,,2 = -K)saw(^i,2) = — J2^^Mn0i,2)- (11) 

^ n=l n 

Fig. 6 shows a comparison between the simple sawtooth function and the numerical result 
obtained using the parameters listed above. This model has also been used to model the 
dynamics of density waves and to calculate the narrow-band noise spectra, with excellent 
agreement with experiment. These results will be reported in some detail in a separate 
pubhcation. 

Eq. (11) is related by charge- flux duality to the current-phase relation of a Josephson 
junction. The periodic behavior of the source-to-drain voltage vs. gate voltage of an STT 
can readily be understood by exploiting the duality with a dc SQUID. When each JJ in 
a dc SQUID has an ideal, sinusoidal current-phase relation and the critical currents Jq are 
identical, the total current is / = Jofsin^Ji + sin 722], where 922 = V^i + 27r$/$o- This yields 
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a total critical current, Ic = 2Iq\ cos(27r$/$o)| , which is a periodic function of the flux $. 
Similarly, the total source-to-drain voltage of an ideal STT, V = V1 + V2, where Vi^2 are given 
by Eq. (11), yields a critical voltage that is a periodic function of gate voltage Vg = Qg/Cg, as 
shown in Fig. 7. Here we note that 62 = 61 + 6g + 6Q, where 6g = 2ni3Qg/Qo, 6q = 27r/3go/<5o; 
Qg is the displacement charge across the gate capacitor, go is the offset charge, and P is the 
screening parameter due to the normal electrons. Fig. 7, although rather simplistic, provides 
at least a preliminary, qualitative interpretation of our experimental results. 

IV. CONCLUSION 

We have carried experiments which, combined with our modeling studies, provide further 
evidence for a quantum mechanism of CDW transport. The use of wave functions, in Eq. 
(7) and (8), is based on the approximation that the 0no-'s are roughly uniform sufficiently 
far from the kink and antikink. An extension of the model would employ wave junction- 
als that incorporated spatial variations of the phases (pnai^)- A soliton domain wall could 
then be viewed as a condensate of a microscopic quantum solitons. The wave functional 
approach to quantum field theory, together with a non perturbative treatment of the func- 
tional Schrodinger equation, has previously been discussed in the context of QCD,0 scalar 



Further refinements might include the addition of disorder to the s-G pinning potential, 
which would enable a description of met ast ability, hysteresis, and related phenomena. Krive 
and Rozhavskii0 have studied the effective Lagrangian for macroscopic quantum tunneling 
of a CDW in a random medium, and have found that disorder significantly enhances the 



included the electrostatic charging energy due to spatial variations of the phase, would 
provide a more realistic description of many aspects of the problem. However, we believe it is 
also important to incorporate the numerous microscopic degrees of freedom when calculating 
the effective action relevant for tunneling. We are thus exploring a generalization of the 
tunneling Hamiltonianil to the case of matrix elements between wave functionals of scalar 
fields representing these degrees of freedom. 

The quantum decay of the false vacuum, accompanied by the nucleation of topological 
defects, is a far-reaching problem, which could impact many areas of physics. For example. 




tunneling rate. A quantum version 



of a modified Fukuyama-Lee-Rice (FLR) model,!! which 
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topological defects, such as flux vortices, play an important role in the cuprates and other 
type-II superconductors.0 Magnetic relaxation rates that depend weakly on temperature up 
to 20 K,il or even decrease with temperature,S suggest that Abrikosov vortices may tunnel 
over a wide temperature range. Moreover, the consistently low IcRn products of cuprate 
Josephson devices suggest that Josephson vortex-antivortex pair creationH may occur when 
the current is much smaller than the "classical" critical current Iq ~ A/Rnc. In cosmology, 
the existence of many matter flelds may facilitate quantum nucleation of a universe even 
when the total action is large, as suggested by Hawking et alS Finally, the extraordinary 
rapidity of flrst-order phase transitions, such as the palpably visible nucleation of ice in 
supercooled water, suggest a possible similarity to the decay of the false vacuum. 
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FIG. 1: (a) CDW phase vs. position, illustrating the production of a soliton-antisoliton domain 
wall pair, (b) Model of a CDW capacitor, showing the nucleated domain walls moving towards 
the contacts. The electric field between the domain walls is reduced by the internal field E*. The 

distances I, Xo, and the crystal thickness are greatly exaggerated for clarity. 

FIG. 2: Potential energy vs. for two different values of 6. Tunneling is prevented by conservation 
of energy when 9 < ir. The small energy barrier for each microscopic degree of freedom within the 
condensate enables tunneling when 9 > n, while the Peierls condensation energy prevents thermal 
hopping. 

FIG. 3: CDW current vs. source-to-drain voltage for several values of gate voltage at 35 K. (The 
shunt current of the normal electrons has been subtracted for clarity.) 

FIG. 4: Vsd vs. Vg for fixed values of total bias current at 35 K. 

FIG. 5: Computed ground state energy E vs. 9, showing the periodic, piecewise parabolic form. 
The apparent divergence for 9 > An is an artifact of the finite number of coefficients used in the 
trial wave function. 

FIG. 6: Voltage as a function of 9, similar to the dual of the Joscphson current-phase relation. 
The sawtooth expansion with 10^ Fourier components (dashed line) is compared to the numerical 
results obtained by minimizing the energy using the trial function discussed in the text. 

FIG. 7: Predicted critical voltage vs. normalized gate voltage 9g for several values of 9o using the 
idealized model discussed in Section III, showing the periodic behavior. 
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TABLE I: Charge-flux duality between density waves and Josephson junctions. 



Density Wave Josephson junction 

Soliton or antisoliton kink w/ charge ibQo Josephson vortex w/ flux ±$0 
Type of threshold Threshold field Et Threshold current It 

Transport characteristic I vs.V, I = 2-KQod(l)/dt Vvs.I, V = 2Tr^od(j)/dt 
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